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g^ ABSTRACT 

< 

We observed the TeV blazar 1ES 1218+304 with the X-ray astronomy satellite 
Suzaku in May 2006. At the beginning of the two-day continuous observation, 
we detected a large flare in which the 5—10 keV flux changed by a factor of ~2 
on a timescale of 5xl0 4 s. During the flare, the increase in the hard X-ray flux 
clearly lagged behind that observed in the soft X-rays, with the maximum lag of 
2.3xl0 4 s observed between the 0.3—1 keV and 5—10 keV bands. Furthermore 
we discovered that the temporal profile of the flare clearly changes with energy, 
being more symmetric at higher energies. From the spectral fitting of multi- 
>■ | wavelength data assuming a one-zone, homogeneous synchrotron self-Compton 

^vq \ model, we obtain B ~ 0.047 G, emission region size R = 3.0 x 10 16 cm for an 

appropriate beaming with a Doppler factor of 5 = 20. This value of B is in good 

agreement with an independent estimate through the model fit to the observed 
"xT ' 

time lag ascribing the energy-dependent variability to differential acceleration 
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timescale of relativistic electrons provided that the gyro-factor £ is 10 . 



Subject headings: BL Lacerate objects: individual (1ES 1218+304) - radiation 



mechanisms: non-thermal - X-rays: galaxies 
c3 ' 



institute of Space and Astronautical Science/JAXA, Sagamihara, Kanagawa 229-8510, Japan 

2 Department of Physics, Tokyo Institute of Technology, Meguro-ku, Tokyo 152-8551, Japan 

3 Stanford Linear Accelerator Center and Kavli Institute for Particle Astrophysics and Cosmology, Stan- 
ford University, Stanford, CA 94305, USA 

4 Universitat Wiirzburg, Am Hubland, D-97074 Wiirzburg, Germany 

5 Landessternwarte, Universitat Heidelberg, Konigstuhl, 69117 Heidelberg, Germany 



-2- 



Introduction 



Blazars are a sub-category of Active Galactic Nuclei where a relativistic jet pointing close 
to our line of sight produces Doppler-boosted emission (e.g., Urry & Padovani 1995; Ulrich, 
Maraschi & Urry 1997). Generally, their overall spectra have two pronounced continuum 
components: one peaking between IR and X-rays is produced by the synchrotron radiation 
of relativistic electrons, and the other in the 7-ray regime, presumably due to the inverse 
Compton (IC) emission by the same electrons. In some cases, 7-ray emission is seen to 
extend to the TeV range; the X-ray and GeV/TeV 7-ray bands correspond to the highest 
energy ends (E max ) of the synchrotron/IC emission (e.g., Inoue & Takahara 1996; Kirk, 
Rieger & Mastichiadis 1998). At these ends, variability is expected to be most pronounced, 
and in fact, such large flux variations are observed, on a timescale of hours to days (e.g., 
Kataoka et al. 2001; Tanihata et al. 2001) or even shorter (minutes scale; Aharonian et 
al. 2007; Albert et al. 2007). Using ASCA data, Takahashi et al. (1996) argued the soft 
X-ray (< 1 keV) variation of Mrk 421, observed to lag behind that of the hard X-rays (> 2 
keV) by ~ 4 ks, may well be ascribed to the energy dependence of the synchrotron cooling 
timescale. More recently, Kataoka et al. (2000) interpreted an observed soft-lag and spectral 
evolution of PKS 2155-304 by a newly developed time-dependent synchrotron self-Compton 
(SSC) model. 

The above paradigm of "soft-lag" was questioned, however, in several aspects. First, 
intensive X-ray monitoring of blazars has revealed not only soft lags but in some cases hard 
lags (Takahashi et al. 2000) which may be a manifestation of another process, e.g., energy 
dependent acceleration. Second, Edelson et al. (2001) voiced concerns about reliability of 
measurement of lags that are smaller than the orbital periods (~ 6 ks) of low Earth orbit 
satellites. This was refuted by Tanihata et al. (2001) and Zhang et al. (2004) who showed 
that, although periodic gaps introduce larger uncertainties than evenly sampled data, lags on 
hour-scale cannot be the result of periodic gaps. A time resolved cross correlation analysis of 
uninterrupted Mrk 421 data obtained by XMM- Newton revealed lags of both signs, changing 
on time scales of up to a few 10 3 s (Brinkmann et al. 2005). Hence the situation is very 
complex and still under debate. 

In this letter we present new results from the May 2006 Suzaku observation of 1ES 1218+304 
conducted as part of a multi-wavelength campaign with KVA-Swi/t-MAGIC. 1ES 1218+304 
is categorized as a high-frequency peaked BL Lac object, at a redshift z = 0.182 (Veron- 
Cetty & Veron 2003). It was discovered as a TeV emitter by MAGIC at energies > 100 GeV 
(Albert et al. 2006) and subsequently confirmed by VERITAS (Fortin 2007). While the 
detailed multiband analysis is ongoing, we focus in this letter on a remarkable X-ray flare 
observed with Suzaku. We present temporal and spectral features in § 3, in § 4 we discuss 
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a physical origin of temporal variability. 



2. Observation and Data Reduction 

1ES 1218+304 was observed with Suzaku (Mitsuda et al. 2007) during 2006 May 20-21 
UT, yielding a net exposure time of 79.9 ks. Suzaku carries four sets of X-ray telescopes 
(Serlemitsos et al. 2007) each with a focal-plane X-ray CCD camera (XIS, X-ray Imaging 
Spectrometer; Koyama et al. 2007) that is sensitive over the 0.3-12 keV band, together with 
a non-imaging Hard X-ray Detector (HXD; Takahashi et al. 2007; Kokubun et al. 2007), 
which covers the 10-600 keV energy band with Si PIN photo-diodes and GSO scintillation 
detectors. 1ES 1218+304 was focused on the nominal center position of the HXD detector. 

For the XIS, we analyzed the screened data, reduced via Suzaku software version 2.0. 
The screening was based on the following criteria: (1) only ASCA-grade 0,2,3,4,6 events 
were accumulated, while hot and flickering pixels were removed using the CLEANSIS script, 
(2) the time interval after the passage of South Atlantic Anomaly is greater than 500 s, (3) 
the object is at least 5° and 20° above the rim of the Earth (ELV) during night and day, 
respectively. In addition, we also select the data with a cutoff rigidity (COR) larger than 
6 GV. After this screening, the net exposure for good time intervals is 69.4 ks. The XIS 
events were extracted from a circular region with a radius of 4.2' centered on the source peak, 
whereas the background was accumulated in an annulus with inner and outer radii of 5.4' 
and 7.3', respectively. We checked that the use of different source and background regions 
did not affect the analysis results. The response and auxiliary files are produced using the 
analysis tools xisrmfgen and xissimarfgen developed by the Suzaku team, which are 
included in the software package HEAsoft version 6.4. 

The HXD/PIN data (version 2.0) were processed with basically the same screening crite- 
ria as those for the XIS, except that ELV > 5° through night and day and COR > 8 GV. The 
HXD/PIN instrumental background spectra were provided by the HXD team for each obser- 
vation (Kokubun et al. 2007; Fukazawa et al. 2006). Both the source and background spectra 
were made with identical good time intervals and the exposure was corrected for detector 
deadtime of 6.0%. We used the response files version ae_hxd_pinhxdnom2_20070914.rsp, 
provided by the HXD team. 
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3. Analysis and Results 

Figure [1] shows the averaged light curves of the four XISs in the six X-ray energy bands. 
Although we could see variations of count rates at some level using HXD/PIN data, it was not 
significant within uncertainties of photon statistics. Thus in the following, we concentrate 
on the temporal variability of the XIS data only, below 10 keV. The temporal variation of 
the hardness ratio (HR) is also shown in the bottom panel of Figure [D It indicates that the 
variability in the soft and hard X-ray bands are not well synchronized. 

To quantify the different shape of the flare with energy dependent time-lags, we fitted 
the light curves with a function given by Norris (1996) after a slight modification of adding 
a constant offset Co to mimic the observed light curves: 

I(t) = C + Ci x exp[-(|i - t pcak /(j r |) fc ] (for t < £ pcak ), 
= C + Ci x exp[-(|t - £pcak/^d|)1 (for t > £ pcak ), 

where £ peak is the time of the flare's maximum intensity Ci, k is a measure of pulse sharpness, 
cr r and <7d are the rise and decay time constants. If the light curve is symmetric in time, 
<7 r and o"d are expected to be equal. All the light curves were binned at 2880 s (a half of 
the orbital period of Suzaku) for fitting. The results of the fittings are given in Table [TJ 
In summary, the observed flare shows the following characteristics: (1) The flare shape is 
asymmetric in time (cr r /o"d < 1) especially in the lower energy band (but note cr T /a^ ~ 1 
for 5—10 keV light curve). (2) The flare amplitude defined as {C\ + C )/C becomes larger 
as the photon energy increases (the 5—10 keV flux changed by a factor of ~2). (3) The 
rise-time of the flare is almost constant ~5x 10 4 s below 2 keV, while it becomes gradually 
longer at higher energy bands. 

Next, we try to evaluate lags of temporal variations in various energy bands. Taking into 
account a wide variety of the flare shape measured at different energies, we estimated lags 
by just comparing the peak-time of the flare rather than using other temporal techniques, 
such as the discrete correlation function (DCF; Edelson & Krolik 1989) or the modified mean 
deviation method (MMD; Hufnagel & Bregman 1992). \j We compared the peak-time in 
five lower energy bands to that determined in the 5—10 keV band. Apparently, the hard 
X-ray (5-10 keV) peak lagged behind that in the soft X-ray (0.3-1 keV) by (2.3±0.7)xl0 4 
s. Importantly, this is much larger than the orbital period of Suzaku and less affected by 
artifacts proposed in Edelson et al. (2001). 



1 Since the DCF quantifies the degree of similarity or correlation between two time series as a function of 
the time-lag, it is not suitable to evaluate "energy-dependent" profiles, as observed in 1ES 1218+304. 
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Fig. 1. — Light curves of 1ES 1218+304 observed with Suzaku XISs in 2006 May. The energy 
bands are 0.3-1, 1-1.5, 1.5-2, 2-3, 3-5 and 5-10 keV (from the upper panel), respectively. The 
bottom panel shows the HR of count rates, defined as (5-10 keV)/(0.3-l keV). The dotted 
line is the constant offset Cq. The dashed line is the characteristic variability time scale of 



a flare t v 



5xl0 4 s. 



The time averaged four XISs and HXD/PIN background subtracted spectra were fitted 
using XSPEC ver. 11.3.2, including data within the energy band 0.6—50 keV. The background 
of HXD/PIN includes both the instrumental (non X-ray) background and the contribution 
from the cosmic X-ray background (CXB; Gruber et al. 1999). Here the form of the CXB was 
taken as 9.0 x 10~ 9 (-E/3 keV) -029 exp(— E/40 keV) erg cm -2 s _1 sr _1 keV -1 and the observed 
spectrum was simulated assuming the PIN detector response to isotropic diffuse emission. 
We first fitted with a single power- law model with Galactic absorption iV H = 1.78 x 10 20 
cm -2 (Costamante et al. 2001). We obtained the best fit photon index T = 2.14 ± 0.01, 
but this model did not represent the spectrum well yielding a reduced x 2 OI 1-23 for 1967 
dof. We also tried to fit with a broken power-law model with Galactic absorption. The 
photon index below the break energy E^ vk (Ti) is 2.04 ± 0.01 while the index above E^^ is 
2.17 ± 0.01, where E hrk (T 2 ) is 1.42 ± 0.05 keV. The flux over 2-10 keV is ~ 2.0 x 10 -11 erg 
cm -2 s -1 . This model gives a better fit with a reduced x 2 of 1.14 for 1959 dof compared to the 
single power-law model, but x 2 is still not acceptable. Considering the spectral variability, 
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Table 1: Tempo ral profiles during the flare of 1ES 1218+304 in various X-ra y energy bands. 
E (keV) t P eak (10 4 s) a r /a d k (d + C )/C ' 



0.3-1 


4.4±0.1 


0.34+0.14 


1.2+0.4 


1.2±0.1 


1-1.5 


4.7+0.3 


0.31+0.08 


1.7+0.5 


1.2±0.1 


1.5-2 


4.7+0.3 


0.31+0.11 


1.4±0.5 


1.2±0.1 


2-3 


5.1+0.3 


0.33+0.10 


1.4±0.4 


1.3±0.1 


3-5 


6.1+0.4 


0.67+0.12 


2.7+0.4 


1.4+0.1 


5-10 


6.7+0.7 


0.84±0.17 


2.8+0.6 


1.6±0.1 



we analyzed the spectrum every 5760 s. The power-law indices vary from 2.05 ± 0.01 to 
2.22 ± 0.01 during the flare, and each segment can be fitted well with a single power-law 
model or broken power-law model with % 2 /dof ranging from 0.94 to 1.09. 

Figure [2] shows the spectral energy distribution (SED) of 1ES 1218+304 with currently 
available datasets. The TeV data are obtained from Albert et a!. (2006), and are corrected 
for the absorption due to the IR Extragalactic background light (EBL; see Figj2]); other data 
are from the NED database. Note that the TeV analysis of multi-wavelength campaign data 
is still ongoing and the combined datasets will be investigated in forthcoming paper (Stefan 
et al. in prep). As expected from the curved X-ray spectrum with photon index T around 2 
and E cnt < 10 keV, the synchrotron emission peaks just around the Suzaku bandpass. 

In order to specify the SED of 1ES 1218+304, we applied a one-zone homogeneous 
SSC model developed in Kataoka et al. (1999). Noting that the characteristic variability 
time scale of the flare is t var ~ 5xl0 4 s, which is most probably determined by the light 
travel time across the source emitting region (see discussion in § 4), we obtain R = ct vav 5 = 
3.0xl0 16 cm for a moderate beaming factor of 5 = 20 (e.g., Kataoka et al. 1999; 2000 for 
self-consistent determination of physical parameters in TeV blazars). With this parameter 
set, the SED of 1ES 1218+304 is fitted with B = 0.047 G, s = 1.7, 7 min =l, 7 br k = 8.0 x 10 3 
and 7 max = 8.0 x 10 5 . We also note that the energy densities of electrons and fields are 
u e = 8.3 x 10~ 3 erg/cm 3 and ub = 8.8 x 10~ 5 erg/cm 3 , respectively. Thus the jet in 
1ES 1218+304 is particle dominated, and the ratio u e /us ~ 100 is well within the range of 
typical TeV blazars. 



4. Discussion 

In our observation we detected a large flare during which the hard X-ray variation lagged 
behind that in the soft X-rays, ~ 2.3 x 10 4 s. This is completely opposite to a well-known 
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Fig. 2. — Overall SED of 1ES 1218+304. Filled circles show X-ray data (Suzaku; this work). 
For TeV data (Albert et al. 2006; filled diamonds), we adopt the correction for the IR EBL 
by Primack et al. 2001 (open diamonds); see also Primack et al. (2005) as well as the 
discussion based on the recent HESS detection of distant TeV blazars (Costamante 2007). 
The other plots are from the NED database. The dashed line is a prediction by a one-zone 
SSC model (Kataoka et al. 2000). 



behavior in which the spectra harden during the phases of rising flux, and soften during the 
phases of falling flux, as has been obtained from the past observations. In the theoretical 
context, however, "hard lag" is actually expected especially in the X-ray variability of TeV 
blazars, but has never been observed so clearly before. It has been suggested that a hard-lag 
is observable only at energies closer to the maximum electron energy, 7 max (Kirk, Rieger & 
Mastichiadis 1998), where the acceleration time is almost comparable to the cooling time 
scale of radiating electrons: £ acc (7max) - £ C ooi(7max)- 

It is convenient to express £ acc and £ coo i in terms of the observed photon energy E (in 
units of keV). Noting that the typical synchrotron emission frequency, averaged over pitch 
angles, of an electron with energy 7771c 2 is given by v ~ 3.7xl0 6 i?7 2 Hz, we obtain; 

t acc (£) = 9.65xl0- 2 (l + z) 3 /^ J B- 3 / 2 r 3 / 2 E 1 / 2 s, 
t cool (E) = 3.04xl0 +3 (l + ^) 1 / 2 5- 3 / 2 r 1 / 2 J E;- 1 / 2 s, 

where z is the redshift, B is the magnetic field strength, £ is the "gyro-factor" which can 
be identified with the ratio of energy in an ordered magnetic field to that in a turbulent 
magnetic field (£ = 1 for the Bohm limit; see, e.g., Inoue & Takahara 1996), and 5 is the 
beaming factor. Note that for lower energy photons, t& cc (E) is always shorter than t coo \(E) 



because higher energy electrons need longer time to be accelerated (£ acc (7) oc 7) but cool 
rapidly (i C ooi(7) oc 7 -1 )- This energy dependence of acceleration/cooling time-scales may 
qualitatively explain the observed characteristics of the X-ray light curves of 1ES 1218+304. 
It is thus interesting to consider a simple toy model in which the rise time of the flare 
is primarily controlled by the acceleration time of the electrons corresponding to observed 
photon energies, while the fall time of the flare is due to the synchrotron cooling time scale. 
In this model, the amount of "hard-lag", r hard , is simply due to the difference of t acc , and 
independent of the energy dependence of t coo \: 

'hard ''accl v -^- / hi/ l acc\-*- J \aw ) 

~ 9.65 x 1(T 2 (1 + zfl\B-*/H-*l\Elf - £>£) s, 



where E\ ow and E^ are the lower and higher X-ray photon energies to which the time-lag is 
observed. Here we took ^low/hi to be the logarithmic mean energy in the observation energy 
bandpass. The result of the model fit to the observed rh ar d is shown in Figure [3] (left). 

Assuming a beaming factor 5 = 20 from multiband spectral fitting (see § 3), the best 
fit parameter of the magnetic field B can be written as ~ 0.049^5 G, where £5 is the "gyro- 
factor" in units of 10 5 . Thus, in order to have the B field required in the acceleration region 
consistent with that derived from the SED fitting, we infer £ ~ 10 5 . Such high value of £ 
is in fact consistent with that inferred by Inoue & Takahara (1996) for other blazars. With 
these parameters, the maximum synchrotron radiation energy -E max , corresponding to 7 max , 
is expected to be ~ 5.3 keV. Hence, the above toy model qualitatively well represents the 
observed spectral/temporal features of 1ES 1218+304, in particular: (1) the synchrotron 
component peaks around the Suzaku XIS energy band in the multiband spectrum (Figure 
[2]) and (2) the observed light curve is symmetric in shape when measured at the high energy 
band, while being "asymmetric" (i.e., fall time longer than the rise time) at the lower energy 
band. Figure [3] (right) compares the energy dependence of observed and modeled flare 
shapes, defined as the ratio of rise and decay time-scales, cr r /cr d . The dashed line shows 
the model prediction from a r /a d ~ t acc /t CO oi = (E/E max ) l ^/(E/E max )-^ 2 ~ E/5.3 keV. 
Although the general trend is well reproduced, Figure [3] indicates that the observed rise 
time may have a bit longer time scale than expected from the model. The most natural 
interpretation for this is the smoothing of rapid variability by the source light crossing time 
scale t crs (e.g., Chiaberge & Ghisellini 1999; Kataoka et al. 2000). Hence if the acceleration 
time scale is shorter than the source crossing time, we expect t crs to smooth out t acc . The 
dash-dotted line in Fig. [3] (right) shows the ratio of the time scales of t crs /t coo i- we can 
see that t crs is longer than t acc below ~ 2 keV, but comparable or shorter above ~ 2 keV. 
As a result, for 1ES 1218+304 it seems reasonable that the rise time of the flare is primary 
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determined by the acceleration time of the electrons, while the fall time of the flare is due 
to the synchrotron cooling time scale. Only the rise time of the flare in the lower energy 
bands are dominated by the light crossing time £ crs , but further studies using more data are 
necessary to confirm this model. 

We are grateful to Dr. L. Costamante for kindly providing us with the numerical 
values of the EBL effects on TeV gamma-ray spectra, and Dr. M. Sikora for many fruitful 
discussions. This work was supported, in part, by a Department of Energy contract to SLAC 
no. DE-AC3-76SF00515. 
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Fig. 3. — Left: Time lag of photons of various X-ray energy bands vs 5-10 keV band photons. 
The solid line corresponds to a fit with r hard = 9Mx I0~ 2 (l + zf^B- 3 / 2 5- 3 / 2 (7. l^-E^) 
where 5 is fixed to 20.0. Right: Energy dependence of the pulse shape, defined as the 
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